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Perturbative expansion of Chern-Simons theory 

Justin Sawon 

An overview of the perturbative expansion of the Chern-Simons path integral is 
given. The main goal is to describe how trivalent graphs appear: as they already 
occur in the perturbative expansion of an analogous finite-dimensional integral, we 
discuss this case in detail. 

81T18; 57M27, 58J28, 81T13 



1 Introduction 

While a standard technique used by physicists, many mathematicians are not familiar 
with the mechanisms and justifications behind perturbatively expanding Feynman 
path integrals. In this article we will describe how trivalent graphs arise when one 
expands a path integral whose Lagrangian contains a cubic term. We'll start with a 
finite-dimensional integral, and then indicate the necessary adjustments required to 
generalize to the infinite-dimensional integral arising in Chern-Simons theory. 

This article is for the proceedings of the BIRS workshop "The interaction of finite 
type and Gromov-Witten invariants". The purpose of the workshop was to investigate 
'large- N duality', which relates perturbative Chern-Simons theory with gauge group 
SU(A0 to certain open and closed string theories. One of the more mysterious aspects 
of this correspondence, to this author at least, is that in string theory trivalent graphs 
(more precisely, fat graphs) have some geometric realization, whereas in perturbative 
Chern-Simons theory they enter purely as part of the technique of calculation. As we 
shall see, trivalent graphs arise merely as combinatorial objects which index terms in 
the series expansion. 

The exact solution to Chem-Simons theory was described by Witten [25] in the late 
1980s. The perturbative expansion of the path integral for flat M 3 was then considered 
by Guadagnini, Martellini, and Mintchev [10, 11, 12, 13], and for general three- 
manifolds by Axelrod and Singer [3, 4]. Subsequently much attention was focused on 
understanding the perturbative series for knots and links in M 3 (see B ar-Natan [5 , 7] , Bott 
and Taubes [8], Altschuler and Freidel [1], and Thurston [24]). More on perturbative 
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expansion of path integrals may be found in Witten's lectures [27]. For a slightly 
different application, to the perturbative expansion of Hermitian matrix integrals, see 
Mulase [16], where Feynman diagram techniques are nicely explained. 

While this article was being prepared, a preprint by Michael Polyak [19] appeared, 
covering very similar material. All the same, the author decided to complete this article 
for inclusion in the present volume. The interested reader may wish to consult Polyak's 
preprint, which in many cases provides additional details to what is covered here. 
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2 Stationary-phase approximation 

Fix a compact, connected, simply-connected gauge group G (for example, SU(n)). 
These conditions imply that tt\(G) = ~ki(G) = 1 and ttt,{G) = Z, a fact that we shall 
use shortly. Let M be a three-manifold without boundary. The field in Chern-Simons 
theory is a connection V on a principal G -bundle over M. The conditions on G imply 
that the bundle can be trivialized and hence we can write the connection as V = d + A 
where A G Q l (M) g is a Lie algebra valued one-form. A change of trivialization is 
given by a gauge transformation, which is a smooth function g : M — > G. Under this 
transformation A changes to g Ag + g~ l dg. 

Definition 2.1 The Chern-Simons action is given by 



Remark Let us explain the terms in the integrand. Write A = ( Ej where £ are 
one-forms and {£/} is a basis for q. On forms, A denotes the usual wedge product. On 
the Lie algebra part, the first term Tr(EjEj) is an invariant bilinear form Bjj normalized 
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so that B(h a , h a ) = 2 , where h a is a coroot corresponding to a long root a . The second 
term looks like 



where cfj are the structure constants of g . We usually write CjjBlk simply as cjjk ■ For 
a matrix group such as SU(ra) , these operations can be performed simply as they appear: 
by multiplying and taking traces of matrices. 

Observe that cjjk gives a three-form wong defined by 



where X, Y, Z € Q and B is the invariant bilinear form as above. This can be made 
into a left-invariant three-form on G, and f y u> = Sir 2 where 7 is a generator of 
H 3 (G,Z) = Z. 

The action S gives a functional on the space of connections A. Let us calculate its 
derivative. 

Lemma 2.2 The curvature F of A can be regarded as a one- form on A, and the 
derivative of S equals ^F. 

Proof Let 5 A describe a variation of the connection A. We can regard 5 A as an 
element of TaA, and then the curvature F of A defines a one-form by 



Tr(A A A A A) = Tt(EiEjE k )C'a ( j A ( k 

= |Tr([£ / ,£ 7 ]^)C / AC 7 AC 



■K 



u(X, Y,Z) := B([X,Y],Z). 




Moreover, 



S(A + <5A) - S(A) = 




M 



Tr((5A A dA + A A d5A + 25A A A A A) 




where we have used integration by parts. 



□ 
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Remark The gauge group Q acts on A, and F is an invariant one-form. On the other 
hand, a calculation shows that 



where deg(g) is the degree of g as a map from M to G, that is, g takes the fundamental 
class of M to deg(g)7, where 7 is the generator of H?,(G, Z) mentioned earlier. Note 
that the degree gives an isomorphism tvq{Q) = Z. 

Since S is well-defined on A/Q up to a multiple of 2ir, the integrand in the next 
definition is invariant under gauge transformations. 

Definition 2.3 The partition function of Chem-Simons theory is given by the path 
integral 



J A 

where the level k is a positive integer. 

Remark The measure on the infinite-dimensional space A of connections is not 
well-defined, so the path integral is really a heuristic tool. It suggests that there should 
be a topological invariant Z^(M) of M. Moreover, the path integral suggests various 
ways of studying this invariant, for instance via a topological quantum field theory or 
via a perturbative series. 

In the perturbative approach k should be thought of as being proportional to the inverse 
of Planck's constant h~ l . In the classical limit K — ► 0, so we should expand in k~ l . In 
the k — > 00 limit the path integral localizes around stationary points: this should be 
understood as follows. Assume that A is not a stationary point of the action S. The 
integrand exp(ikS(A)) lies on the unit circle, and if we vary A slightly the value of the 
integrand rotates around the circle; in the large-& limit the values that the integrand 
takes are distributed evenly around the unit circle and so cancel out in the integral. 

By Lemma 2.2 the derivative of S vanishes at flat connections, so these are the stationary 
points. To expand around a flat connection a, we let A = a + j3 where j3 G Q l (M) (g) g 
should be regarded as being the perturbation from a. A calculation shows that 



where d a is the connection on one-forms induced by the connection a . The contribution 
to Zji{M) from the flat connection a is therefore 



S(g- i Ag + g~ l dg)-S(A) 



2vrdeg(g) 





(1) 
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where the first factor exp(ikS(a)) in the integrand is a constant and the second factor 
involves both a quadratic and a cubic term in (5. Note that the linear term in [3 vanishes 
because we are expanding around a stationary point. 

This process of localization reduces the path integral to an integral over the moduli 
space of flat connections. If the flat connections are isolated then Z/,(M) is simply a 
sum of integrals like the one above. However, the flat connections are never isolated 
because of the action of the gauge group - we will return to this point later. Note also 
that (1) is still an infinite-dimensional integral; our first step will be to calculate the 
integral in the simpler finite-dimensional case. 

3 A finite-dimensional model 

We will consider the following finite-dimensional integral 



where Q and V are respectively quadratic and cubic forms on Mr (linear in all entries 
and totally symmetric), and k is a positive real constant. One could drop the symmetry 
assumption on V, which would require some additional labeling of the half-edges at 
trivalent vertices of the graphs arising later, but in this article we will only consider 
symmetric V, which is why we have included the factor of 1/3! . 

The integral is not absolutely convergent, but it can be defined in some formal sense 
via analytic continuation. To begin with, the Gaussian integral on R can be evaluated 





for k real and positive. Analytic continuation in k yields 




In dimension greater than one, we can diagonalize the quadratic form to get 



A" 




for Q positive definite, and 
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for Q indefinite (but still non-degenerate). We will refer to all of these integrals, which 
don't involve cubic terms, as Gaussians. 

Returning to Z*, we first deal with the dependence on k by a change of variable 
x i — y := \fkx (we immediately rewrite the dummy variable x 1 as x in our formula) 

Z k = Jfc"? / d N xe i{kQ^)+i^vi W )) 

a OO .jyj 

= IT? / d N xe il 2 Q(x ' x) Y" — — w V(x,x,x) m . 
Jm* ~ m\6 m k2 

We have used the power series for the exponential function involving the cubic term, 
but notice that the odd terms in the series give an integrand which is an odd function on 
R N , so those terms do not contribute to the integral. On the other hand, the 2m th term 
gives the coefficient of k~ 2 m as 

Jw (2m)!6 2m 
To calculate this we need to understand how to integrate the Gaussian multiplied by a 
polynomial. 

Lemma 3.1 Let J(x) = JjX> be a linear functional on M. N (we are using the Einstein 
summation convention) and define 

Z(J) := / d N xe i ^ Q(x ' x)+J{x) ) 

so that Z(0) is just the Gaussian integral. Then 

Z(J) = Z(0)e- L 2 Q ' i{J ' J) 
where Q~ l (J,J) = Q> k JjJk and Q k is the matrix of the inverse of Q . 

Proof This is a simple exercise in completing the square. □ 

The linear term / is known as a source, and the main trick is to differentiate Z(7) with 
respect to J and then evaluate at / = . We will illustrate this through some examples. 

Example First observe that, using the definition of Z(7) , 

£/Wi««V- (-^)zwU„. 

Now using the expression for Z(J) in Lemma 3.1, the right-hand side becomes 

-iZ(V)(-i&%e-&- l V>y j=0 . 
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This vanishes when we evaluate at J = , as it should because it is the integral over M. N 
of an odd function. 

Next observe that to integrate the product of the Gaussian and a degree two polynomial, 
we should differentiate Z(J) twice 

jL'^-(-i)(^K 

Once again, we use Lemma 3.1 to evaluate this as 

(-ifZ(O)(-i0 k e-^- l(J > J) + J k (- ■ .))| /=0 = 2(0)0* 

The precise expression (• • • ) is irrelevant as it vanishes when we evaluate at / = 0. 

In general, to integrate the product of the Gaussian and a degree d polynomial, we 
should differentiate Z(J) d times. 

Now let us apply this to the m = 1 term occurring in Zt 

Writing V(x,x,x) in coordinates as Vjkix^x k x l or V a b c x a x b x c , we see that the integrand 
is the product of a Gaussian and a degree six polynomial, and hence the integral equals 



-1 



2!6 2 



7=0 

V jkl V ahc i 3 Z(O)(0 k Q ,a Q bc + Q ik Q lb Q ac + ■■■) 



The last factor on the right consists of the sum over all ways of dividing the set of 
indices {j, k, I, a, b, c} into pairs, and hence contains jt ( 2 2 2 ) = ^ terms. Any other 
terms which appear when we differentiate Z(/) will vanish when we evaluate at / = . 

Since we are using the summation convention, j,k, etc, are dummy indices. Up to a 
factor of (^)Z(O), we get nine terms of the form V j idV a bcQ! k Q la Q hc and six of the 
form Vj k iV ah cQ>' a Q kh Q lc ■ These terms can be encoded using the trivalent graphs 

and T? = 



More precisely, to each trivalent graph V we associate a weight W(T) given by labeling 
vertices with V and edges with Q 1 , and then contracting indices. For instance on 
Ti we first label the legs (ends of the edges) with indices. Then we label the vertices 
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with Vju and V a bc , and the edges with Qj k , Q la , and Q hc , in such a way that the indices 
agree with those on the legs. 




Qla 




Qik 



Then we apply the summation convention and arrive at the weight 
(2) W(Ti) := V jkl V ahc Q k Q la Q hc . 

In this construction V and Q~ 1 are known as the vertex and propagator, respectively. 
Summing up, the m = 1 term becomes 

1 ^jZ(Q)(9W(Y l ) + 6W(Y 2 )) = iZ(Q)(\W(Y x ) + ^W(Y2)) 

= iZ ^V^\ W{Tl)+ \^\ W{T2) 

where AutTi and Autl^ are the automorphism groups of the graphs. We think of a 
trivalent graph as a collection of half-edges partitioned into two element subsets (edges) 
and three element subsets (vertices). Then in this context, an automorphism of the graph 
means a permutation of the half-edges which preserves both the partition into edges and 
the partition into vertices (this differs from the usual definition). For example, Aut Y\ 
is generated by the involutions 

(D-Q 0±0, and CMD, 

whereas Aut T2 is the semi-direct product of the group generated by the involution 



and the symmetric group S3 permuting the three edges. 

Using the same arguments, we deduce that the 2m th term looks like 

A_^_ V ^...V fl6c i 3 -Z(0)(fi*...fi fc + ---) 

with 2m copies of V, and the last factor being the sum over all ways of dividing a set 
of 3 x 2m indices into pairs, thus a sum of j^y ^^™ 2) terms - Each trivalent graph Y 
with 2m vertices gives us a weight W(Y) as before, which we claim occurs precisely 
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^Autr" ti mes m tne sum above. To understand why, it helps to separate the trivalent 
vertices from the edges of the graph. 



2m < 



V 



V 



<2 1 



> 3m 



D or 

Then Y is given by gluing together the legs in some specific way. Note that if the legs 
are labeled with indices, then gluing corresponds to contracting indices. Now observe 
that the semi-direct product of groups S2 m x (Si) 2 " 1 acts on the trivalent vertices. Under 
this action, we still get Y but the dummy indices change, unless the element of the 
group induces an automorphism of Y . Thus the orbit has size ( ^utr| ' corresponding 
to the number of times W(Y) occurs in the sum above. 

Thus the 2m th term is 



; m z(0) £ 



l 



(m+l) — loop 
graphs 



AutT 



W(Y) 



(trivalent graphs with 2m vertices are usually referred to as (m + l)-loop graphs), and 

1 



Z k = Y,i m Z(0)k- N >- m 



m=0 



(m+l) — loop 
graphs 



AutT 



W(Y). 



This sum over graphs includes disconnected graphs. However, the multiplicativity 
of W(Y) and the simple behaviour of AutT under disjoint union readily leads to the 
following formula: 

x ' ffi=Q connected 



connected 
(m+l) — loop 
graphs 



4 Degenerate quadratic forms 

Before returning to the infinite-dimensional case, we need to consider how to modify 
our formulae in the case that the quadratic form is degenerate. Recall that the Gaussian 
integral for a non-degenerate form Q is given by 



/ 



^ ' xe ik?Q(x,x) 



l„ r , / 2 ^\ f e ^'(signG)/4 



k 



det£| 
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Now suppose that a compact group K acts on M. N preserving Q. We can write 

f d N xe ikk 2 QM = [ d(ie ikl 2 Q(x ' x) Vo\O x 
ii* Jr n /k 

where dp is the induced measure on the space of orbits R N /K, and VolO x is the 
volume of the orbit through x G M. N . More precisely, the integral over R N /K can be 
thought of as an integral over a slice for the -action, that is, a submanifold of R N 
intersecting the generic orbit in exactly one point. The action induces a map 

p x : K -► O x C R N 

taking g G K to g(x) £ O x . Denote the derivative of this map by 

B:=d Px : t • T x O x C T X R N = R N 

where t is the Lie algebra of K. Then 

VolO* = / drf 

drj\detB\ 
= Vol^|detB*B|^ 

where drj is the Haar measure on K, dr( is the induced measure on O x , and B* is the 
adjoint of B. Since p x is ^-equivariant, detB is constant on K, justifying the last 
equality. We usually refer to |detB*B|2 as the Jacobian determinant. Summarizing, 
we have 

(3) f d N xe ikL 2 Q{x ' x) = Yo\K [ dfie ikl 2 Q(x ' x) \detB*B\ 1 2. 

Jr n Jr n /k 

Example The usual S l -action on M 2 preserves the standard quadratic form x 2 + y 2 . 
In this case, maps 9 to (— Oy, 9x) and the adjoint B* x ^ maps (X, Y) to —Xy + Yx, 
so 

\detB*B\ k 2 =( X 2 +y 2 ) 1/2 = r 

and thus 

/ dxdye~( x +y ^ = 2tt / <ire _r r = 7r. 

In the example above, Q is still non-degenerate. However, if the action of the group 
K is affine, then this forces the form Q to be degenerate in the direction of the orbits. 
The author has been unable to find such an example with K compact, though in 
infinite-dimensions the action of the gauge group on the space of connections is affine. 
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Now when K is non-compact, a A"-invariant function will be non-integrable. However, 
the integral on the right-hand side of (3) can still be finite, and we will use this as a 
definition of the left-hand side (we will also ignore the infinite Vol A' factor in this case). 

We assume that Q induces a non-degenerate form on the space of orbits. Thus we 
can reduce the original ill-defined integral involving a degenerate quadratic form to 
one over the space of orbits, with the inclusion of a factor proportional to the Jacobian 
determinant. This additional factor can be dealt with via the Faddeev-Popov method, 
which we will describe shortly, but first we will compute the one-loop term directly in 
the infinite-dimensional case. 

5 The one-loop term 

In this section we will calculate the one-loop term in perturbative Chern-Simons theory 
following Atiyah [2]. Recall that the contribution to the Chern-Simons path integral 
coming from connections A = a + in a neighbourhood of the flat connection a is 



In this section we will ignore the cubic term in (3 and just focus on the infinite- 
dimensional Gaussian integral coming from the quadratic term, which is known as the 
one-loop term. In other words, we will evaluate 



First observe that after choosing a Riemannian metric on M, we can define a pairing on 
n\M) (8) by 



where * is the Hodge star operator acting on the one-form part of 02 (which of course 
depends on the choice of metric). Recall that for X and Y G su(n), — Tr(XY) is a 
positive-definite invariant bilinear form, which explains why we have included a minus 
sign in the above pairing. The quadratic form Q is then given by the self-adjoint 
operator — * d a on Q l (M) tg) a, as 



This is completely analogous to the finite-dimensional case, where the quadratic form is 
given by a matrix combined with the standard pairing on R . 
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Now the operator Q = — * d a is degenerate on fr(M) © q, and the analogue of 
B : t — > T x O x C is the map coming from infinitesimal gauge transformations 

In other words, Q is degenerate on the image d a £l (we use fi 7 ' as an abbreviation for 
SV(M) © 0). If we assume that the flat connection a is irreducible and isolated, then 
the twisted de Rham complex 

fi°(Af) ^(M) © 2 (M) © fl\M) © 

is exact. Specifically, a flat connection gives a representation of 7Ti(M) which is 
irreducible if and only if H°(M, 0) vanishes, and H l (M, g) is the space of first-order 
deformations of the connection, which vanishes if the connection is isolated. The 
remaining cohomology groups H 2 and H 3 are dual to H 1 and H° respectively. 

According to the previous section, our Gaussian integral is defined as 

I dPe ik ^Q'^\detB*B\-2. 

The quadratic form Q' is given by the reduction of Q on the quotient Q l /d a Q°; note 
that it is non-degenerate by exactness of the de Rham complex. 

Define an operator P on 

n° © n 1 = ft © d a n° © (n 1 /^ ) 

by (d a , d^, — * d a ) where d* a := — * d a * is the adjoint of d a on one-forms. Writing 
this as a matrix (with respect to the three factors) 

/OB* \ 
P = \ B 
V a I 

we see that 

|detP| = |detS*B||det(2'|. 

Moreover, 

B*B = A° (the Laplacian on Q°(M) © 0) 

P 2 = A°ffiA^ (the Laplacian onfi°(M)©0©17 1 (^) ®0)- 

Note that although Q is degenerate on Q l , we can extend it to a non-degenerate 
quadratic form P on the larger space 0° © O 1 . (This partly explains the addition of the 
ghost field c £ f2° in the Faddeev-Popov method - see the next section.) 
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The Jacobian determinant | det B*B\ 2 = (det A°)2 is independent of (3, while when 
we integrate the Gaussian we pick up a factor of | det Q \ 2 . This combination can be 
rewritten 



detfi*fi|2 



det B*B 



det Q'\2 



det PI 2 



(detA°)* 
(detAL)* 



For the determinants of the Laplacians, we need to calculate the product of their 
eigenvalues, which can be done using zeta-function regularization. This gives the 
one-loop term, up to a phase factor and an overall constant. 

Bearing in mind that the Chern-Simons partition function is supposed to represent a 
topological invariant of M, one might wonder about the appearance of a Riemannian 
metric in our calculations above. While it is true that the Laplacians (and their 
determinants) depend on the choice of metric, Ray and Singer showed that the ratio 



is independent of this choice [20]. They showed this by explicitly calculating the 
variation under a change of metric, and this quantity became known as the Ray-Singer 
analytic torsion. It was later shown that it equals the Reidemeister torsion, which 
is defined combinatorially. In particular, our factor above is the square root of the 
Ray-Singer torsion, and is therefore a topological invariant. 

6 The Faddeev-Popov method 

There is another way of dealing with the Jacobian determinant, which is necessary to 
extend our results to the higher order terms in the perturbative expansion. This is the 
method used by Axelrod and Singer [3] and Bar-Natan [5, 7], and involves introducing 
additional fields and integrating over them in the path integral, as we now explain. 

To begin with, we want to integrate over the space of orbits of the gauge group acting 
on Q l (M) <S> q. One way to achieve this is to find a function F whose level sets meet 
each orbit transversely in a single point, and then integrate over F~ 1 (0) . Of course, we 
must also include the Jacobian determinant, which represents the volume of the orbit. 



(detA°p 
(detA^ 
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Equivalently, we could include a 5-function term 5(F(f3)) in the path integral. Up to a 
factor, this can be written as 

5(F(P)) oc J d^ m '^ 

by using the Fourier representation of the J-function. The term F{j3)cf) can then be 
added to the Lagrangian, and becomes an additional field to integrate over. In the 
Chern-Simons case, F comes from the gauge-fixing condition d^/3 = 0, and <p should 
be a g-valued function. 

The Jacobian determinant can also be written in terms of additional fields. For this we 
must introduce ghost fields c and c, which lie in f2°(M) (g) g but are Fermionic, that is, 
anti-commutative (whereas f3 and <fi are Bosonic). Then, up to a factor, 

|detfi*fi|3 = (detA°)^ oc J dc dc e ik ^ A °<* c \ 

which is the analogue of the Gaussian in the theory of anti-commutative integration. 
So our new path integral looks like 

J d(id^dcdce ikS{f3 ^ c) 

where after the addition of gauge-fixing and ghost terms the action is given by 

5(/3,0,c,c) = — / Tr(/3 A d a (3 + |/3 A A (3 + d* /3 A *<j> + c A *A° c) . 
47r Jm 

One might be concerned about the appearance of metric dependent terms in the new 
Lagrangian, such as the *-operator and Laplacian. To show that the path integral is 
independent of the choice of metric we introduce the BRST operator. In this subsection 
only, Q will denote the BRST operator, which should not be confused with the quadratic 
form Q defined earlier. It is an odd (interchanges Bosonic and Fermionic fields) 
derivation on the space of functionals in /3, (p, c, and c defined by 

QP = -{d a + (3)c 

Q4> = 

Qc = <t> 

Qc = |[c,c]. 

When we deform the metric, only the gauge-fixing and ghost terms vary. However, one 
can show that this variation is Q-exact. Moreover, Q is like a vector field with zero 
divergence, so the integral of a Q-exact term must vanish. See Axelrod and Singer [3] 
or Bar-Natan [7] for the details. 
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The introduction of the BRST operator can perhaps be understood as an algebraic 
method of forming a quotient, though as far as the author knows, this idea has not yet 
been mathematically formalized 1 . In our case, we are interested in the quotient of the 
space of connections by gauge transformations. Note that Q extends the infinitesimal 
gauge action on (3 to the other fields 4>, c, and c. It is straightforward to check that 
Q 2 = ; essentially what we are doing above is taking Q-cohomology. Only the metric 
invariant part 

^- [ Tr(/3 A d a f3 + 1 13 A p A 0) 
47r Jm 

survives to give an element of H°(0, which is our algebraic representation of the space 
of functionals on the quotient. The BRST approach to gauge theory is described in 
Witten's third lecture in [26]; Chapter 7 of Bar-Natan's thesis [5] is another useful 
reference. 

It should also be noted that a more natural way to view the above path integral is in its 
superspace formulation. However, since we wish to write the terms of the perturbative 
expansion as configuration space integrals we must stick to the language of differential 
forms. Both of these approaches are described by Axelrod and Singer [3]. 



7 Higher order terms 

We will now give a rough outline of how one constructs the higher order terms in 
the perturbative expansion of the Chern-Simons partition function. These terms will 
be indexed by trivalent graphs, and formally the series will look the same as the 
finite-dimensional case, once we have properly interpreted the quadratic and cubic 
forms on S7 ! (M) (g) q. Actually, the series has only been fully computed in the case that 
a is the trivial connection, in which case the twisted de Rham complex decouples 

Vt\M) ®g A CT(M) ®g A VL Z {M) <g> q Vl\M) ® g 

where d is the usual exterior derivative on forms on M, and Id is the identity on g. We 
will assume a is trivial for the remainder of this article; of course if M is Euclidean 
space then the only flat connection is the trivial connection. 

We saw that the quadratic form Q is given by the self-adjoint operator — * d a , when 
combined with the pairing 

{Pi, Pi) :=~ I Tr(A A*/3 2 ) 
27 r Jm 

'Note added in proof: for a homo logical algebraic interpretation due to Jim Stasheff see 
[21,22, 23]. 
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on O (M) (8) g. When a is the trivial connection this decouples into — * d (combined 
with the induced inner product) on Q 1 (M) and the invariant quadratic form fi on g. 
Recall that to calculate weights on trivalent graphs we needed the inverse QT 1 . As far 
as g is concerned, this is just the inverse B~ l of the invariant quadratic form. 

On Vt\M) we need to calculate the inverse of the elliptic operator — * d. Strictly 
speaking, this inverse does not exist as — * d has non-trivial kernel. However, the kernel 
is the image of the infinitesimal gauge action, and is taken care of by the Faddeev-Popov 
method described in the last section. After adding gauge-fixing terms, we find that the 
new Lagrangian has a quadratic term which looks like 



where L := — (*d + d*). Moreover, L 2 is the Laplacian (on functions and one-forms), 
whose inverse is given in terms of a Green's function. For example, in Euclidean space 
IR 3 the solution to Aw = w is given by 



The Green's function is 



and the inverse of L itself is given by 

T „ e^(x; - yi)dxj A dx k 

LoG = — . ^ 

4tt\x — y\ 3 

(using Einstein summation), meaning that this is the kernel of an integral operator 
representing L 1 . 

The cubic term V is much simpler. Recall that if we write (3 = ('Ej where C 7 are 
one-forms and {£/} is a basis for g, then 



where cuk comes from the structure constants of g. So the cubic term decouples into 
simply wedging the three one-forms together (and then integrating over M) on Q 1 (M) 
and cuk on g. 

The Faddeev-Popov method also introduces ghost fields, so there is an additional 
propagator and interaction coming from ghosts. These are usually drawn with dashed 
lines in our Feynman diagrams, whereas the propagator and interaction described above 
are drawn with solid lines. See Bar-Natan [5, 7] for more details. 



ik\(P + <l>,Ltj3 + <!>)) 




Tr(/3 A A 0) = \cukC 1 A C 7 A C 



K 
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Let r be a trivalent graph. Since both the quadratic and cubic terms decouple, the 
weight W(T) will become a product 

(4) W(T) = I r (M)xb s (T) 

where Ir(M) does not depend on the Lie algebra g and b 5 (F) does not depend on the 
three-manifold M. The fact that a is the trivial connection is crucial for this decoupling, 
as otherwise the quadratic form would not decouple. 

Remark For each finite-dimensional Lie algebra g , b s is a weight system on trivalent 
graphs. It is the same as the weight system arising in the theory of Vassiliev invariants [6]; 
for example 

b s (T 2 ) := c UK c ABC B IA B JB B KC . 

Notice that the cubic terms on Q}(M) and g are both skew-symmetric. Consequently 
we really need to introduce the notion of an orientation on T (an equivalence class of 
cyclic orderings of the legs at each trivalent vertex), so that Ir(M) and b g (F) make 
sense. The skew-symmetry cancels out in W(F): if we reverse the orientation on T, 
both Ir(M) and b g (T) will change sign but W(T) will remain the same. 

One consequence of this skew-symmetry of vertices is that many weights automatically 
vanish. For example, write Ti and Ti for the graph 




with its two orientations. Then Ti and Ti are actually isomorphic as oriented graphs, 
so 

b B (Ti) = b B (Ti) = -6 fl (ri) 

must vanish. 



8 Configuration space integrals 

We will describe in more detail the terms Ir(M) arising from (4). In our finite- 
dimensional model we summed over indices when evaluating W(T), as in (2), and this 
still works fine for the Lie algebra factor b g (T) since g is finite-dimensional. However, 
on the infinite-dimensional part summation is replaced by integration: heuristically, 
one can regard the ^-functions {5(x — y)} y eM as a basis for functions on M, so that 
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summation over elements in the basis becomes J yeM (and something similar happens 
for one-forms). The term Ir 2 (M), for example, therefore looks like 

dx dyK(x,y) A K(x,y) AK(x,y). 



M JM 



In this formula, K(x, y) is the Green's function for L on M. We can only write this 
out explicitly when M is Euclidean space, since we have an explicit formula for the 
Green's function on M 3 . 

For a graph T with m vertices, we get an integral over M m . Now the Green's function 
K(x, y) will have a singularity at x = y, so the integrand only makes sense on the 
configuration space 



However, by analyzing the severity of the singularity Axelrod and Singer [3] showed that 
the integrals on M m are actually finite. One approach here is to extend the integrand to 
the compactification of (M) due to Fulton and MacPherson [9], which is a manifold 
with boundary and corners, and then show that this extension is finite. 

It is important to note that for a fixed trivalent graph r , Ir is not a topological invariant 
of three-manifolds: it varies as we deform the metric on M. The variation can be 
understood as an integral over the boundary of the Fulton-MacPherson compactification 
of C°(M). However, when we combine the terms I? in the sum 



the variations should cancel out. In fact there are two kinds of variation, and variations 
of the first kind cancel due to a special property of the weights b & (T), namely that they 
satisfy the IHX relations. Variations of the second kind are known as "anomalous", 
and are more difficult to deal with. If we choose a framing of M, then we can correct 
the anomaly with a counter-term and thereby construct a genuine invariant of M (see 
Sections 5 and 6 of Axelrod and Singer [3]). 

9 Invariants of knots and links 

Finally let us say a few words about the generalization to knot and link invariants. Given 
an oriented knot K in M, we can include an additional term in the Chern-Simons path 
integral given by the monodromy of the connection A around /C . The monodromy 



C°(M) := {(* (1) , . . . ,x (m) ) G M m | x {i) + x® for i + j}. 



(5) 




graphs 



Qeometry & Topology Monographs 8 (2006) 



Perturbative expansion of Chern-Simons theory 



163 



would take values in g, so we evaluate its trace in a representation V of G. This gives 
the Wilson loop 

W V (X) :=Tr y Pexp / A. 

Jkcm 

In this expression, the path ordered exponential Pexp can be represented by a sum of 
iterated integrals 

oo « « 

(6) Y, I ' ' ' / A (1C(h)) ® • • • ® A(£(fc)) 

* =0 0<fl<f 2 <-<f*:<l 

where we have used /C : S 1 — > M to denote a parametrization of the knot. We then 
represent the g-valued A(/C(f;)) in End(V), multiply them, and take the trace. The 
expectation value of the Wilson loop 

Z k {M;K)= [ VAe ikS(A) W v (IC) 
J A 

should be an invariant of the knot. For links one includes several Wilson loops, one for 
each component of the link. 

The expectation value Zk(M; K) admits a perturbative expansion just like Zjt(M), and 
once again each term splits into the product of a weight system (coming from the Lie 
algebra g and representation V) and a configuration space integral. The Feynman 
graphs consist of unitrivalent graphs, whose univalent vertices lie on a circle S 1 (or 
collection of / circles in the case of an /-component link). The interactions at the 
univalent vertices come from the terms A(/C(f,-)) in (6). The configuration space looks 
like 

C%(M; JQ := {(x (1 > . . . , x (a )x {a+l ] x {a+b) ) G M" x K, b \ x w + x ij) for i + j} 

and the integral can be replaced by an integral over a suitable compactification of 
C° ab (M;K.) (see Thurston [24]). 

A great deal of work has gone into understanding these configuration space integrals 
when M is Euclidean space, since we then have an explicit formula for the Green's 
function of L _1 (see Bar-Natan [5, 7], Bott and Taubes [8], Altschuler and Freidel [1], 
and Thurston [24]). The simplest example comes from the graph given by a single 
chord connecting two circles, in the case of a two component link C (it looks like Ti , 
but with the circles corresponding to the two components of the link). Then 

Cj {11} (M;£) = Ci x C 2 = S l x S l 

Qeometry & Topology Monographs 8 (2006) 



164 



Justin Sawon 



is already compact (we write {1,1} to indicate that there is a point on each component 
of the link). The integral is given by 



and it computes the degree of the map from C\ x C 2 = S l x S l to S 2 which takes (x, y) 
to j^E^r- It i s known as the Gauss linking number of C\ and £2- 

The Gauss linking number is a genuine invariant, but for a general configuration space 
integral there are two kinds of variation when we isotope the knot or link. It is easy to 
show that, when combined in a sum like (5), the variations of the first kind cancel out. 
However, the second kind of variation is more difficult to deal with, and once again 
leads to an "anomaly". 

For example, the first-order term 7e(IR 3 ; /C) for a knot /C is known as the writhe. It 
corresponds to the graph given by a circle with a single chord (it looks like T2 , but with 
the circle corresponding to the knot), and can be expressed as a Gauss integral similar to 
the one above. This integral varies as we isotope the knot, and can take all real values, 
so there is an anomaly. However, if the knot is framed, that is, a trivialization of its 
normal bundle is given, then we get a genuine invariant by adding the total torsion r 
(see Pohl [17]). In fact, up to a factor, 



is the linking number of the knot with its parallel (the knot obtained by using the framing 
to displace K, slightly). 

It was shown independently by Altschuler and Freidel [1] and by Thurston [24] that, up 
to a correction by the anomaly, the perturbative series gives an invariant of knots and 
links in M 3 . This extended earlier work by Bar-Natan [5, 7] and Bott and Taubes [8], 
who had investigated the series up to the degree two term. Moreover, the perturbative 
series was shown to be a universal Vassiliev invariant (c.f. Bar-Natan [6]). It is expected 
that it should agree with the universal Vassiliev invariant of Kontsevich [14], which 
would follow from the vanishing of the anomaly in degree greater than one. This is still 
an open conjecture, though recently Poirier [18] showed that the anomaly does indeed 
vanish in degrees two to six (see also Lescop [15]). 
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